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Abstract 

We study Higgs boson pair production processes at future hadron and lepton colliders including 
the photon collision option in several new physics models; i.e., the two-Higgs-doublet model, the 
scalar leptoquark model, the sequential fourth generation fermion model and the vectorlike quark 
model. Cross sections for these processes can deviate significantly from the standard model predic- 
tions due to the onedoop correction to the triple Higgs boson coupling constant. For the one- loop 
induced processes such as gg — > hh and 77 — > hh, where h is the (lightest) Higgs boson and g and 7 
respectively represent a gluon and a photon, the cross sections can also be affected by new physics 
particles via additional one-loop diagrams. In the two-Higgs-doublet model and scalar leptoquark 
models, cross sections of e + e~ —> hhZ and 77 — > hh can be enhanced due to the nondecoupling 
effect in the one-loop corrections to the triple Higgs boson coupling constant. In the sequential 
fourth generation fermion model, the cross section for gg — > hh becomes very large because of the 
loop effect of the fermions. In the vectorlike quark model, effects are small because the theory has 
decoupling property. Measurements of the Higgs boson pair production processes can be useful to 
explore new physics through the determination of the Higgs potential. 
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I. INTRODUCTION 



The standard model (SM) for particle physics has experienced a great success in describing 
the experimental data of high energy physics below the energy range of a hundred GeV, 
but its portion for electroweak symmetry breaking, the Higgs sector, remains unknown. 
Experimental confirmation of the Higgs boson is one of the most important issues in the 
high energy physics. The direct search results at the LEP experiment have constrained the 
mass (m h ) of the Higgs boson as m h > 114.4 GeV [lj in the SM with one Higgs doublet, and 
the global analysis of precision measurements for electroweak observables has indicated that 
rrih is smaller than 157 GeV at the 95% confidence level [H, E|. In addition, the combined 
data from the CDF and DO collaborations at the Fermilab Tevatron have excluded the region 
of 162 GeV < m h < 166 GeV Q 1 . The CERN Large Hadron Collider (LHC) has already 
started its operation, and it will soon be ready for hunting the Higgs boson. We expect that 
the Higgs boson will be discovered in coming several years. 

Once the Higgs boson is found at the Tevatron or the LHC, its property such as the mass, 
the decay width, production cross sections and the decay branching ratios will be thoroughly 
measured as accurately as possible in order to confirm whether it is really the particle 
responsible for spontaneous breaking of the electroweak symmetry. The Higgs mechanism 
will be tested by determining the coupling constants of the Higgs boson to the weak gauge 
bosons. The measurement of the Yukawa coupling constants will clarify the mass generation 
mechanism of quarks and charged leptons. However, in order to understand the physics 
behind the electroweak symmetry breaking, the Higgs potential must be reconstructed by 
measuring the triple Higgs boson coupling constant (the hhh coupling constant). 

On the other hand, from the theoretical view point, it would be expected that the SM is 
replaced by a more fundamental theory at the TeV scale. One way to see the new dynamics 
is to measure effective vertices of the SM fields and to compare them to the theoretical 
calculation of radiative corrections. The effect can be significant in the electroweak theory 
especially when the mass of a new particle comes mainly from the vacuum expectation value 
(VEV) of the Higgs field like chiral fermions. In such a case the decoupling theorem js| does 
not necessarily hold, so that the new physics effects do not decouple and are significant. It is 
well known that the systematic study of nondecoupling parameters in radiative corrections 
to the gauge boson two point functions has played an important role to constrain new physics 
models by using the precision data of electroweak observables at the LEP and the Stanford 
Linear Collider (SLC)Q. 

Such nondecoupling effects of new physics particles can also be very significant in the 
radiative corrections to the hhh coupling constant Q. Quartic powerlike contributions 
of the mass of a new particle can appear in the one-loop correction to the hhh coupling 
constant, which can give a large deviation from the SM prediction. For example, in the two- 
Higgs-doublet model (THDM), the hhh coupling constant of the lightest (SM-like) Higgs 



Recently, the bound on the Higgs boson mass from Tevatron experiments has been updated 
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boson can be deviated by (9(100)% due to nondecoupling effects of extra scalar bosons in 
radiative corrections without contradiction with perturbative unitarity It is known that 
such a large deviation in the hhh coupling constant from the SM value can be a common 
feature of the Higgs sector with the strong first order electroweak phase transition 0-0 
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Therefore, the 



which is required for a successful scenario of electroweak baryogenesis 
measurement of the hhh coupling constant at collider experiments can be an important 
probe into such a cosmological scenario. The one-loop contributions to the hhh coup ling 
constant can also be very large in the model with sequential fourth generation fermions 14 



and a class of extended supersymmetric SMs 15 



At the LHC, the measurement of the hhh coupling constant would be challenging. In the 
SM, the cross section of double Higgs boson production from gluon fusion, gg — > hh [16 



can be O(10) fb for rrih = 120-160 GeV with the collision energy of y/s = 14 TeV, 
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and vector boson fusion 



while those of double- Higg s-strahlung qq — > V* — > hhV\i 
qq — > V*V*qq — > hhqq [16J, [l8|, |20J are much smaller. The double Higgs boson production 
mechanism from gluon fusion has been studied in Ref. 2l( with the h — > WW^*' decay 
mode. They conclude that the luminosity of 3000 fb" 1 is required to measure the hhh 
coupling constant at the 20-30% level [21]. For a light Higgs boson (rrih < 130 GeV), the 



main decay mode is h — > bb which cannot be useful due to huge QCD backgrounds, so that 
the hhh coupling constant cannot be accurately measured at the LHC. 

At the International Linear Collider (ILC), the accuracy for measuring the hhh coupling 
constant would be better than that at the LHC depending on the mass of the Higgs boson. 
A Higgs boson pair can be produced in the double-Higgs-strahlung process e + e~ — > hhZ 22 
and the W fusion mechanism e + e~ — > hhvv 23[. At the first stage of the ILC where e + e~ 
energy is 500 GeV the hhh coupling can be measured via the double-Higgs-strahlung process 
for m/j < 140 GeV 24N27|. The evaluation of the statistical sensitivity for the hhh coupling 



constant is about 20% accuracy [24j, |25j. Detailed simulation studies for this process are 



ongoing, which shows that the sensitivity may be lower 28[. The photon linear collider 



(PLC ) op tion may also be useful to explore the hhh c oup ling constant for 120 < rrih < 200 



pup i 

3. 



GeV [29|, |30|. A simulation study is also in progress [3JJ. At the second stage of the ILC 
(t/s = 1 TeV) or the Compact Linear Collider (CLIC) where the collision energy would be 
at a multi-TeV scale, the double Higgs boson pair production from W boson fusion becomes 
important because the cross section is larger due to the t-channel enhancement 25|, |27|. The 



statistical sensitivity for the hhh coupling constant is less than 10% 25 



In this paper, we study how the hhh coupling constant affects cross sections for the double 
Higgs boson production processes gg — > hh, e + e~ — > hhZ, e + e~ — > hhvv and 77 — > hh in 
various new physics models such as the THDM, models with scalar leptoquarks, the model 
with the chiral fourth generation fermions and the model with vectorlike quarks. Cross 
sections for these Higgs boson pair production processes are evaluated, and can deviate 
significantly from the SM predictions due to the deviation in the one-loop corrected hhh 
coupling constant. In these processes, the effect of the deviation in the hhh coupling constant 
mainly appears in the interference of the diagram with the hhh coupling constant and the 
other diagrams. Thus, the sign of the deviation can be important. Also, in the one- loop 
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induced processes such as gg — > hh and 77 — > hh, cross sections can depend on new physics 
particles in additional one-loop diagrams. In the THDM and scalar leptoquark models, 
cross sections for e + e~ — > hhZ and 77 — > hh can be enhanced due to the nondecoupling 
effect on the hhh coupling constant through the extra scalar loops. In the chiral fourth 
generation model, cross sections of double Higgs boson production processes can become 
significantly large, because new particles mediate in the leading order loop diagram as well 
as the nondecoupling effect on the hhh coupling constant. In models with vectorlike quarks, 
the effect on the cross sections are small because of the decoupling nature of the theory. 
By measuring these double Higgs boson production processes at different future collider 
experiments, we would be able to test properties of new physics particles in the loop, which 
helps identify the new physics model. 

In Sec. II, effects of the hhh coupling constant in Higgs boson pair production processes 
gg hh at LHC, e + e~ — > hhZ and e + e~ —> hhvv at ILC and CLIC, and 77 — > hh at their 
photon collider options are discussed. Model dependent analyses for these processes are 
given in Sec. Ill for the THDM, the scalar leptoquark models, the chiral fourth generation 
model, and the vectorlike quarks. In Sec. IV, summary and discussions are given. 



II. THE HIGGS BOSON PAIR PRODUCTION PROCESSES AT COLLIDERS 



In this section, we discuss Hig gs b oson pair production processes gg — > hh 16Nl8|. 



e + e- hhZ 22 



e + e — > hhvv 23 



29| in various new physics models. 



and 77 — > hh 

These processes contain the hhh coupling constant so that they can be used to determine 
the hhh coupling constant at future collider experiments. The effective ggh and 77/1 vertices 
would be precisely measured in the single Higgs boson production processes as gg —> h at 
hadron colliders 32] and 77 — )■ /i resonance production at the PLC 33J, which will be used 
to extract the hhh coupling constant from the one-loop induced processes such as gg — > hh 
and 77 — > hh. In this section, before going to the discussion on the calculation for the cross 
sections in each model, we first consider the results in the SM with a constant shift of the 
hhh coupling constant by a factor of (1 + Ak); 



A 



hhh 



A hhh\ l 



Ak) 



(1) 



where A|^ = ~3m\/v at the tree level 2 with v (~ 246 GeV) being the VEV and 
being the mass of the Higgs boson h. This constant shift can be realized when there is 



2 At the one-loop order, the effective hhh vertex function have been evaluated as [3] 



^hhh(Si m hy m h) 



3ml 



1 - 



iTr 2 v' 2 mf i 



l + O 



(2) 



where N c (—3) is the color factor. The full expression of the vertex function Y^h{p\,P2,pl) is also given 
in Appendix A for completeness. In numerical analysis, we include the SM one-loop correction to the hhh 
coupling constant. 
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FIG. 1: The double Higgs boson production process gg 
collider. 



hh via gluon fusion at the hadron 



the dimension six operator in the Higgs potential 29|, |34j . Quantum corrections to the hhh 
coupling constant due to the bosonic loop can also provide the constant shift of the hhh 
coupling constant approximately 

At the LHC, the largest cross section of the Higgs boson pair production comes from the 
gluon fusion mechanism 16Ml8l|. Feynman diagrams for gg — > hh are depicted in FIG. [H 
The triangular loop diagrams contain information of the hhh coupling constant. The parton 
level cross sections are calculated at the leading order as [13] 



d(gg -> hh) 



dt 



1 1 1 



2al 



2 2 8 2 2! 16tts 2 (4tt) 2 



KhhV 



F A + F L 



mt 



+ \G L 



(3) 



where is the loop function for the triangular diagram, while Fq and Gu are those for box 
diagrams which, respectively, correspond to the invariant amplitudes for same and opposite 
polarizations of incoming gluons 21|: see Appendix B. The invariant mass distribution can 



be obtained by multiplying the gluon-gluon luminosity function as 



where M-, 



hh 



da(gg — > hh) 
dMh~ h 

\fi, t = s/s, and 



a(gg -> hh)—^. 

s dr 



(4) 



dL 



gg 



dr 



1 dx 



x 



fg(x,fM F = M hh )f g (r/x,fi F = M hh ), 



(5) 



where f g (x,fi F ) is the parton distribution function of gluons. In our numerical calculation, 
the CTEQ6L parton distribution function is used 35|. The loop integrals are evaluated by 
a package; Loop Tools 
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It is well known that this process receives large QCD corrections 3 . Although the NLO 
calculation is very important in evaluating this process, throughout this paper we totally 



3 The next-to-leading order (NLO) QCD corrections to this process have been computed in the heavy 
top-quark mass limit in Ref. 37(, which give an over all factor K ~ 1.9 (if- factor) for fj, F = Mhh- The 

The correction mainly comes from the 
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smaller value of K ~ 1.65 for [ip = was suggested by Ref. 
initial state radiation of gluons. It is known that this kind of approximation works well in the single Higgs 
production via the gluon fusion mechanism, where the NLO cross section is evaluated by the leading order 
gg — > h cross section for a finite top-quark mass with the if -factor in the large m t limit. The running of 



strong coupling constant can also change the cross section by 25-50% [21] 
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FIG. 2: The invariant mass distribution of the cross section of gg — > hh process at the LHC with 
y/s = 14 TeV for m h = 120 GeV (left) and m h = 160 GeV (right). The solid, dotted, dashed, 
long-dashed and dot-dashed curved lines denote the SM prediction, the SM with the positive 100% 
correction to the hhh coupling constant, that with the +20% correction, that with the —20% 
correction, and that with the —100% correction, respectively. 



neglect NLO QCD corrections in our calculations of the cross section in various new physics 
models. The QCD corrections in each new physics model are currently unknown so that the 
computation of these corrections is beyond the scope of this paper. 

In FIG. [21 we show the invariant mass distributions of the cross section of gg — > hh 
process with the deviation of the hhh coupling constant for = 120 GeV (left) and 
for mh = 160 GeV (right). Throughout this paper we take the top-quark mass to be 
171.2 GeV. These solid, dotted/dashed and long-dashed/dot-dashed curves represent the 
SM prediction including the SM one-loop effect on the hhh coupling constant, that with 
constructive deviations An = +1.0 and +0.2 and that with destructive deviations Ak = —1.0 
and —0.2, respectively. The total cross section is about 20 (10) fb for mh = 120 (160) GeV 
in the SM. Only for An = +1.0, a small peak comes from the large hhh coupling constant 
through the triangular diagram in the near threshold region. The peaks can be found around 
Mhh ~ 400 GeV, which are caused by the interference effect of the triangular and the box 
diagrams. Since these two contributions are destructive to each other, the positive (negative) 
variation of the hhh coupling constant makes the cross sections small (large) in this process. 
This means that in the gg — > hh process the sensitivity is getting better for the negative 
contribution to the hhh coupling constant and vice versa. If we have additional colored 
particles in the new physics model, this situation could be changed. 

At an electron-positron linear collider, the hhh coupling constant will be measured by the 
double-Higgs-strahlung 22j and the Higgs boson pair production via the W boson fusion 
mechanism 23j. Feynman diagrams for these processes are shown in FIG. [21 The e + e~ — > 
hhZ process may be a promising channel at the ILC to measure the hhh coupling constant 
for light Higgs bosons because of the simple kinematical structure. Since relatively larger 
collision energy is required for three body final states of hhZ, the s-channel nature of the 
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FIG. 3: The double Higgs boson production at the e + e collider. The double-Higgs-strahlung 
process e + e~ — > hhZ and the vector boson fusion process e + e~ 
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FIG. 4: The cross sections of e + e —> hhZ process at the ILC as a function of collision energy y/s 
for m h = 120 GeV (left) and m h = 160 GeV (right). 



process may decrease the cross section. On the other hand, if we have large enough energy, 
one can control the collision energy to obtain the maximal production rate. In FIG. HI the 
cross sections of the double-Higgs-strahlung are evaluated as a function of e + e~ center of 
mass energy ^/s. The left (right) panel shows the case with the Higgs boson mass to be 
m/j = 120(160) GeV. The curves are presented in the same manner as in FIG. |2j Under the 
variation of the hhh coupling constant, the cross section of the double-Higgs-strahlung has 
the opposite correlation to that of gg — > hh. Therefore, the positive contributions to the 
hhh coupling constant has an advantage to obtain better sensitivities. 

At a high energy lepton collider, the hard photons can be obtained from the Compton 



back scattering method 38]. By using hard photons, Higgs boson pairs can be produced in 



77 — > hh process. Feynman diagrams for this process are shown in FIG. El and the helicity 
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FIG. 5: The double Higgs boson production process 77 — > hh at the photon collider. 



specified cross sections are given by 



a 



AlA2 = ?( 7 A 1 7A 2 ^ /l/l) 



dt 



2 

a EM 



t -_ 2! 16vrs 2 {An) 2 



Xhhh v 



H 



A1A2 



s — m 



2 ""A 
h 



rrAiA2 



(6) 



where H^ Xi and H^ lX ' 2 are the loop functions 29] (see Appendix B). The total cross section is 



calculated by convoluting with the photon luminosity function / 7 (y, x), where x = 4E e u)o/ml 
can be controlled by the frequency Uq of the laser photon, as 



a 



yy—thh 



Vm fVm dy 

dr I — 

V 



1 + CiC 2 



7 £-7 



a 



1 - 



■a 



f 7 (y,x)fJr/y,x), 



(7) 



where £ 7 is the mean helicity of the photon. The photon luminosity spectrum is given by 



D(x) 



D(x 



i-y 

4 8 , 
1 =) ln(l + x 



1-y- 4r(l - r) - 2\ e X 1 rx{2r - 1)(2 - y) 



(8) 



x x' 



+ 2A e A < 



2 x 2(1+ a;) 2 
5 1 1 



(1 + -)ln(l + x, 

v x 1 v ' 2 1 + x 2(l + x) 2 _ 



(9) 



where r 



x(l-y) 



and A e (A 7 ) is the helicity of the incident electron (photon) 38|. The 



maximal energy fraction of photon y m = is fixed by the kinematics of the Compton 
scattering at the photon collider. In FIG. El the full cross sections of e~e~~ — > hh 
are shown as a function of the energy of the e~e~ system. We here choose the same sign 
polarizations for initial photons in order to efficiently extract information of the hhh coupling 
constant. The parameter x is taken to be 4.8, which can be tuned by the frequency of the 
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FIG. 6: The full cross section of e e (t(+)t(+)) — > hh process as a function of yfs ee for = 120 
GeV (left) and m h = 160 GeV (right). 
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FIG. 7: The cross sections of e + e — > hhvv process at the ILC as a function of collision energy yfs 
for m h = 120 GeV (left) and m h = 160 GeV (right). 



laser photon. The curves are given in the same manner as in FIG. [21 The situation is very 
different from gg — > hh at the LHC Energies of initial gluons are widely varied at a hadron 
collider, while back-scattered photons at the PLC have narrow band spectra. Therefore, we 
can tune the effective energy of photons at the PLC to some extent. The relative strength 
of the W boson and the top-quark loop diagrams strongly depends on the collision energy 
and the Higgs boson mass. Only for mh = 120 GeV, the large hhh coupling constant case 
(An = +1.0) shows a peak at the near threshold regime. It is found that the negative 
deviation of the hhh coupling constant makes cross section large for m h = 120 GeV (left), 
while it has an opposite effect on the cross section for = 160 GeV (right). 

If we go to further high energy e + e~ colliders, the second stage of the ILC or the CLIC, 



the Higgs boson pair production via the W boson fusion mechanism becomes important 23 



The cross section increases for higher energy because of the t-channel enhancement of 
W + W~ — > hh subprocess. In FIG. [7J we evaluate the production rate for e + e~ — > hhvv 



by CalcHEP [39j. For both m h = 120 GeV (left) and m h = 160 GeV (right) cases, the 
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cross section simply scales as a function of energy and can be much larger than those of 
e + e~ — > hhZ and 77 — > hh. The An dependence in the cross section of e + e~ — > hhvv is 
opposite to that in e + e~ — > hhZ; i.e., a larger cross section for e + e~ — > hhvv is obtained for 
a larger \Ak\ value with a negative sign. 

III. COMPARISON OF THE HIGGS BOSON PAIR CREATION PROCESSES IN 
DIFFERENT MODELS 



In this section, we study cross sections for the double Higgs boson production processes 



gg hh [16Hl8 



e+e" -> hhZ 22 



29 in four dif- 



e + e~ — > hhvv |23| and 77 — > hh 
ferent models; i.e., the THDM, the model with scalar leptoquarks, that with chiral fourth 
generation quarks and leptons, and that with vectorlike quarks. In these models, one-loop 
contributions to the hhh coupling constant can be nondecoupling even when new particles 
are heavier than the electroweak scale, so that large deviations in the hhh coupling constant 
can affect the cross sections of these double Higgs boson production processes. Unlike the 
analysis with the constant shift with An in the previous section, the energy dependencies 
in the hhh vertex function are also included in our evaluation here. Furthermore, for the 
one-loop induced processes such as gg — > hh and 77 — > hh, the contribution of additional 
one- loop diagrams where the new particles are running in the loop can be significant. 



A. Two-Higgs-doublet model 



The THDM |40| is the simplest extension of the Higgs sector in the SM, which can ap 



pear in various new physics scenarios such as the minimal supersymmetric SM 



41 



the top 



ro lor model El radiative seesaw models for neutrinos QQ and models of *ctroweak 
baryogenesis [44]]. The Higgs scalar doublets interact with other fields purely by the elec- 
troweak force. Therefore, in addition to the change in the cross section due to the quantum 
correction to the hhh coupling constant, the contribution of the charged Higgs boson loop 
can affect the cross section of 77 — > hh, while the loop effect of extra scalar bosons appear 
only through the correction to the hhh coupling constant for gg — > hh, e + e~ — > hhZ and 
e + e~ — > hhvv. 

The potential of the THDM with a softly-broken discrete Z2 symmetry is given by 



V'thdm = ml$l$i + ml$l$ 2 - 
+ A 3 ($ t 1 <l'i)(^$ 2 ) 



(ml§\§ 2 + H.c. 



A 4 ($ T 1 $ 2 )(^<D 1 ) + 



^1 /^.t^. \2 ^ 



(10) 



where 



= 1,2) are scalar isospin doublet fields with the hypercharge of +1/2, which 
transform as $1 — > $1 and $2 —^2 under the Z 2 - Although m\ and A5 are complex in 
general, we here take them to be real assuming the CP invariance. The two Higgs doublet 
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fields can then be parameterized as 



-j=(vi + hi + izi) 



(11) 



There are 8 degrees of freedom in the two complex scalar doublet fields. Three of them 
are absorbed as the longitudinal components of the weak gauge bosons. The remaining 
five convert into the mass eigenstates, two CP even Higgs bosons (h,H), a CP odd Higgs 
boson (A), and a pair of charged Higgs bosons (H^). The eight parameters m\-m\ and 
A]— A5 are replaced by the VEV v, the mixing angle of CP even Higgs bosons a, the ratio 
of VEVs tan/3 = i^/fi, the Higgs boson masses m h ,m H ,m A ,m H± and the soft breaking 
parameter M 2 = m§/ (sin (3 cos (3). The parameters in the Higgs potential can be constrained 
by imposing theoretical assumptions such as perturbative unitarity 45|, |46( and vacuum 
stability 47] . If we require stability of the theory below a given cutoff scale A imposing the 
conditions of vacuum stability and triviality, the Higgs boson parameters are constrained as 



a function of A by the renormalization group equation analysis 48[ 

Here, we consider the "SM-like" case with sin(/3 — a) = 1 where only the lighter CP-even 



Higgs boson h couples to the weak gauge boson as VVh [49(. In this case, all the coupling 
constants of h to the SM particles take the same form as those in the SM at the tree level. 
The hhh coupling constant is also described by the same tree- level formula as in the SM. 
The difference appears at the loop level due to the one-loop contribution of the extra scalars. 
Under the imposed softly-broken discrete symmetry, there can be four types of Yukawa 
Although in general there can be large phenomenological differences 
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interactions 

among the different types of Yukawa interaction, especially in flavor physics, we do not 
specify the type of Yukawa interaction in this paper, because there is no proper difference 
in the discussion here in the SM-like limit where h behaves as if it were the SM Higgs boson 
at the tree level 4 . 

The extra Higgs bosons have been searched at the LEP experiment. The lower mass 
bound for the CP-even Higgs boson is m H > 92.8 GeV, and that for the CP-odd Higgs boson 
is m A > 93.4 GeV in the minimal supersymmetric SM whose Higgs sector is the THDM [lj]. 



The bound for charged Higgs boson mass has also been set as m H± > 79.3 GeV [lj. The 
electroweak precision data from the LEP experiment may indicate that the Higgs sector 
approximately respects the custodial SU(2) symmetry 53j. This symmetry becomes exact 
in the Higgs potential in the limit of m A = m H ± with arbitrary sin(/3 — a) or in the limit 
of m H = m H ± with sin(/3 — a) = 1 (or m h = m H± with cos(/3 — a) = 1). The breaking 
of the custodial symmetry in the two-Higgs-doublet potential gives large contribution to 
the T parameter which is proportional to the mass differences of extra Higgs bosons. In 
order to suppress these contributions, we take their masses to be degenerate in the following 
discussion; i.e., ran = rriA = rn H ±. 



4 The allowed regions of m H ± and tan f3 can receive constraint from the flavor physics data such as b — > S7 
depending on the type of Yukawa interaction 50 , 5^] . 
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FIG. 8: The rates for one-loop contributions from H, A, in the THDM to the hhh coupling 
constant for mh = 120 GeV (left) and for mh = 160 GeV (right). 



The one-loop correction to the hhh coupling constant in the THDM is evaluated as 



pTHDM 

1 hhh 
pSM 
L hhh 



III 



67T 2 v 2 mf l 



M 2 



Hi 



4 



\2'K 2 v 2 m 2 i 



M 2 



m 



H 



4 



127r 2 v 2 m\ 



M 5 



m. 



(12) 



where sin(/3 — a) = 1 is taken. The deviation from the SM results can be very large when 
M 2 ~ 0. The full calculation of the vertex function is shown in Ref. {?} 5 - In FIG. El the 
deviation in the effective hhh coupling constant from the SM value is shown as a function of 
y/§, the energy of h* —> hh in the THDM Q . The mass of the SM-like Higgs boson is taken 
to be m h = 120 GeV (left) and m h = 160 GeV (right). The masses of extra Higgs bosons are 
taken to be m$ = 200 GeV (dotted line), m$ = 300 GeV (dashed line), and m$ = 400 GeV 

m H ±. These effects can be about 120-70% 



(long-dashed line), where m $ = m H = m A 



for rrih = 120-160 GeV with m $ ~ 400 GeV. FIG. |8] shows that the deviation in the hhh 
coupling constant can be approximately described by the analysis with a constant shift by 
the factor of (1 + An). 

In FIG. EE the invariant mass distribution of the differential cross section for gg — > hh at 
the LHC is shown in the THDM 6 . The curves are given in the same manner as in FIG. [HI 
and the SM predictions are also denoted by solid curves for comparison. Higgs bosons do 
not couple to gluons at the tree level, so that the one-loop effect of the extra Higgs bosons 
only appear in the correction to the hhh coupling constant. For larger extra scalar masses 
m$ = 400 GeV, peaks can be found in the near threshold region of a Higgs pair, which come 
from the enhancement of the hhh coupling constant. There are also peaks around M^h ~ 400 



5 In the case other than sin(/3 — a) — 1, the hhh coupling constant can deviate from the SM value at the 
tree level because of the mixing between h and H. The discussion at the one-loop level with including 



such a mixing effect is also given in Rcf. 
6 Cases without sin(/3 — a) = 1 were considered in Ref. 



particles were also discussed in Ref. [55 



54 1 . The s-channel resonance effect of new physics 
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FIG. 9: The invariant mass distribution of gg — > hh process at the LHC with yfs = 14 TeV for 
m h = 120 GeV (left) and m h = 160 GeV (right) in the THDM. 
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FIG. 10: The cross section of e + e — > hhZ process as a function of yfs for mh = 120 GeV (left) 
and m h = 160 GeV (right) in the THDM. 



GeV, which are interference effects between the triangular and the box diagrams. Those 
contributions weaken each other, and hence the enhancement of the hhh coupling constant 
decreases the cross section as in the SM with constant deviation. 

In FIG. [101 we show the cross section of the process e + e~ — > hhZ as a function of the 
collision energy yfs in the THDM 7 . The curves are presented in the same manner as in 
FIG. EO Relatively large nondecoupling effect of the extra scalar bosons can appear in the 
radiative correction to the hhh coupling constant. 

In FIG. ITU the cross sections of the Higgs pair production at the PLC are given for 



the THDM |57H59|. The extra Higgs boson can contribute to the corrections of the hhh 



coupling constant as well as 77 — > hh process. The hhh coupling constant can be probed 
by choosing the collision energy to be near threshold region for relatively heavy extra Higgs 
bosons > 400 GeV. There are threshold enhancement from the box diagrams after 



7 More general types of double Higgs-boson production processes at e + e colliders were studied in Ref. [56 1. 
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FIG. 11: The cross section of 7(+)7(+) — > hh process at the photon collider option at the ILC as 
a function of the e~e~ collision energy for rrih = 120 GeV (left) and mh = 160 GeV (right) in the 
THDM. 
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FIG. 12: The cross sections of e + e — > hhuv process at the ILC as a function of collision energy 
y/s for m h = 120 GeV (left) and m h = 160 GeV (right) in the THDM. 



2m H ±. The details are shown in Ref. 58 



In FIG. CE21 the cross sections for e + e~ — > hhvv with the one-loop corrected hhh coupling 
constant due to extra scalars are shown. As we show in FIG. El the hhh coupling constant 
can deviate from the SM prediction significantly for Mhh ^ 2m$, while for > 2m$ the 
deviation becomes small where Mhh varies from 2rrih to y/s. We find that the large corrections 
in low Mhh region can enhance the cross section by a factor of a few in magnitude. Although 
the positive one-loop correction decreases the cross section, this process is still important 
because the total cross section can be larger than those in other Higgs pair production 
processes. 
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B. Scalar Leptoquarks 



We next consider contributions to cross sections from scalar leptoquarks [60|. Unlike the 
case with the extra Higgs scalar doublet, the scalar leptoquarks are colored fields. They, 
therefore, can affect the ggh and the gghh vertices at the one-loop level, by which the cross 
section of gg — > hh in this model can differ from the SM prediction in addition to the effect 
of the deviation in the hhh vertex. 

We here introduce a complex scalar, leptoquark, 0^ = (3, 1)1/3 or (3, 1)4/3, as an exam- 
ple for such theories, where SU(3)c, SU(2)l and U(l)y quantum numbers are shown. The 
most general scalar potential can be written as 

/ 2\ 2 

Vm(<S>, 0zq) = A h$| 2 - V -\ + Mjj |^| 2 + Xiq |<^| 4 + A' |<^| 2 |$| 2 , (13) 

where $ is the SM Higgs doublet. The mass of leptoquarks is given by tt?^ = Mj^ + ^f-. 

The searches for the leptoquarks have been performed at the collider experiments at LEP, 
HERA and Tevatron. In order to avoid large contributions to lepton flavor violating pro- 
cesses, the leptoquarks are usually assumed to be coupled only with one fermion generation 
in their mass eigenbasis. Under this assumption, the experimental bounds are evaluated 



as ra^ LQ > 256 GeV [6l|, 316 GeV (62[, and 229 GeV [63| for the leptoquarks interacting 



only with the first, second, and third generation, respectively. There are indirect limits for 



masses and their Yukawa couplings through the effective four-fermion interaction [64 . 

In the large mass limit of an SU(2)l singlet scalar leptoquark, the coupling strength of 
the effective ggh vertex is enhanced by the factor of 4/3 as compared to the SM prediction, 
and hence the cross section of the single Higgs boson production gg — > h can be enhanced 
approximately by 16/9. The SM Higgs boson with mass 162-166 GeV has been ruled out by 
analyzing gg — > h — > process at the Tevatron j^j. The exclusion band of the Higgs 

boson mass from Tevatron results can be translated a wider range as 155-185 GeV. 

The one-loop correction to the hhh coupling constant due to the scalar leptoquark can 
be calculated analogous to the charged Higgs boson contribution in the THDM [7| as 

N c mi ( M? n \ 3 
thML ~ 1 + c ^ 1 _ ^ (141 

nf h ~ + 6^ml { ml L J ■ [ } 

The full expression of the one-loop corrected vertex T^ h (s, m\, ml) is given in Appendix A. 
In FIG. [131 we evaluate the relative size of the one- loop contributions to the hhh coupling 
constant from the leptoquarks for = 120 GeV (left) and for = 185 GeV (right). 
Three reference values are taken for leptoquark masses as rn^ LQ = 256 GeV (dotted line), 
m <t>LQ = 300 GeV (dashed line), and m^, 1Q = 400 GeV (long-dashed line). For rrih = 120 
GeV the one-loop correction can be about 150% by a singlet leptoquark with m^, iQ = 400 
GeV, while for = 185 GeV it can be about 60%. These effects are constructive to the SM 
value in the nondecoupling region Mjg ~ 0, which can be significant for heavy leptoquarks 
because of the m\ LQ enhancement. Below the thresholds of 2m^ LQ the quantum effects are 
approximately flat for the function of the off-shell Higgs boson energy. 
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FIG. 13: The rates for one-loop contributions of an SU(2) singlet scalar leptoquark to the hhh 
coupling constant for mh = 120 GeV (left) and for m,h = 185 GeV (right). 
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FIG. 14: The invariant mass distribution of gg — > hh process at the LHC with y/s = 14 TeV for 
rrih = 120 GeV (left) and = 185 GeV (right) in the leptoquark model. 

In FIG. [HI we show the invariant mass distribution of the cross section for gg — > hh 
process at the LHC. For rrih = 120 GeV with relatively larger leptoquark masses m$ > 300 
GeV, peaks can be found in the near threshold region of a Higgs boson pair, which come from 
the enhancement of the hhh coupling constant as well as that of the ggh coupling constant 
due to the new colored particles in the triangular diagram: see FIG. [1] for the reference. The 
threshold enhancement of the on-shell leptoquark pair production can also be seen around 
M hh ~ 2m <j) . 

In FIG. [151 the cross sections for e + e~ —> hhZ are shown in the case with rrth = 120 
GeV as a function of the collision energy y/s for various rn^ LQ . The effects of the scalar 
leptoquarks only appear in the hhh coupling constant. 

In FIG. HEl we show the cross section of the Higgs pair production process at a PLC. 
The effects of the scalar leptoquarks depend on not only their masses but also their electric 
charges. The first peaks around the threshold region come from the modifications of the 
hhh coupling constant and the effective 77/1 vertex due to the scalar leptoquarks. The 
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FIG. 16: The cross section of 7(+)7(+) — > hh process at the photon collider option at the ILC as 
a function of the e~ "e~ collision energy for rrih = 120 GeV (left column) and rrih = 185 GeV (right 
column). The electric charges for scalar leptoquarks are taken as Q = 4/3 (top) and Q = 1/3 
(bottom), respectively. 
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FIG. 17: The cross sections of e + e — > hhuv process at the ILC as a function of collision energy 
yfs for nih = 120 GeV (left) and = 185 GeV (right) in leptoquark models. 

former effect does not depend on the quantum number, while the latter does. The threshold 
enhancements of the box diagram can be found after ^fs^ e ~ 2m^ which are also dependent 
on the electric charges. 

In FIG. [TTJ, the cross section for e + e~ — > hhvv is shown for = 120 GeV (left) and 
rrih = 185 GeV (right) with the hhh coupling corrections of leptoquarks. The cross section 
becomes smaller as in the THDM because of the negative interference. 



C. Chiral fourth generation 

One of the fundamental questions in the SM is the number of generation (family) of quarks 
and leptons. There is no theoretical reason to restrict the fermion families to be three. The 
electroweak precision data also do not exclude completely existence of the sequential fourth 
generation. The fourth generation quarks are colored particles so that they, similarly to 
the case of scalar leptoquarks, affect the ggh and the gghh vertices. In addition, the chiral 
fermion has the nondecoupling property; i.e., the mass is completely proportional to the 
VEV, so that the one-loop correction to the hhh coupling constant can be very large in 
magnitude. 

We here introduce a sequential set of fermions, i.e., Q' = (t' L , b' L ) T , L' = 
(£' L , v' L ) T , t' R , b' R , £' R , i/ R , as the chiral fourth generation model (Ch4). Neutrinos are as- 
sumed to be Dirac particle whose masses are generated by Yukawa interaction. From LEP 
results, a lower bound on the extra charged lepton (£') mass is 100.8 GeV, while for heavy 
neutral lepton (z/) with Dirac nature to be 90.3 GeV. The fourth generation up-type quark 
(t') is rather stringently constrained, m t / > 256 GeV, by Tevatron. This bound is inde- 
pendent of the Cabibbo-Kobayashi-Maskawa (CKM) elements between the first three SM 
fermions and the fourth generation. A similar limit for the heavy down-type quark (b') is 



obtained as my > 128(268) GeV by using the charged (neutral) current decay modes 64 



Further stronger bounds can be found in Ref. 65( by assuming additional assumptions. 
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(S,T) 


m t ' — my = 50 GeV 


m t i — my = 55 GeV 


m t / — my = 60 GeV 


m t , = 256 GeV 


(0.18,0.18) 


(0.18,0.22) 


(0.17,0.26) 


m t i = 300 GeV 


(0.19,0.18) 


(0.18,0.22) 


(0.18,0.26) 


mt , = 400 GeV 


(0.19,0.18) 


(0.19,0.22) 


(0.19,0.26) 



TABLE I: The contributions from the chiral fourth generation fermions to S and T are shown. 
The mass degeneracies m t r = m^ and my = m v > are assumed. 



The contributions from the fourth generation fermions to the oblique electroweak pa- 
rameters can be significant. In the limit of heavy fermions, the S"-parameter is calculated 

as 

2 1 f, m\, ml,' 



S = In -4 - In . (15) 

3vr 6tt V ml ml 1 



It is noted that a complete set of the four generation fermions with degenerate masses is 



excluded at 6a level, if we take into account the constraint only from the S'-parameter 641.166 



It can be relaxed by requiring the preferable mass hierarchy of the fourth generation fermions 



i.e., m t ' > my and my > m v t [I4j. Furthermore, the breaking of isospin symmetry in the 
fourth generation doublet gives substantial positive contribution to the T-parameter, which 
pulls back the model to the allowed region in (S,T) plane. For m^ = 117 GeV, the 95% 



CL upper bounds are given as S < 0.16 and T < 0.21 with a strong correlation 64j, [66 



In Table (TJ the contributions from the fourth generation fermions to oblique parameters are 
listed. In order to reduce the number of parameters, m t > = m^ and my = m v < are taken. We 
can see that the appropriate mass difference (m t i — my = 50-55 GeV) between the fourth 
generation fermions fits the constraint from the electroweak precision data. 

From a theoretical point of view, there are a few constraints such as triviality bounds 
and vacuum stability bounds on the coupling constants 67J]. To avoid the Landau pole for 
the fourth generation Yukawa coupling constant, the mass should be lighter than about 577 
GeV when the cutoff scale of theory is taken to be 2 TeV. Instability of the vacuum gives 
more serious bound on the Higgs boson mass. It may be evaded by introduction of the extra 
Higgs doublet 67J or some other new physics dynamics, and hence we here keep the mass 
of the (SM-like) Higgs boson to be the electroweak scale. 

The additional heavy colored particles enhance the effective ggh coupling approximately 
by a factor of 3 which leads to enhancement of the cross section of go_-^ h by a factor of 
9 at hadron colliders. Tevatron bounds on the SM Higgs boson mass [3j can be translated 
into a wider exclusion band as 125-200 GeV for the Higgs boson mass in the chiral fourth 



generation model 68 



The large one-loop effect on the hhh coupling constant in the SM is generalized straight- 
forwardly to the chiral fourth generation model: 

hhh ~ 1 - \ / (16) 

hhh f'=t',b',e'y h 
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FIG. 18: The rates for one-loop contributions of the chiral fourth generation fermions to the hhh 
coupling constant for nih = 120 GeV (left) and for = 210 GeV (right). The dotted, dashed, 
long-dashed curved lines indicate masses of heavy fermions as m t > = 256, 300, 400 GeV, respectively, 
with the appropriate mass difference nit* — my = 55 GeV. 

The explicit formula of including energy dependence is given in Appendix A. Since rrif, 
enhancements come from extra heavy fermions, we would expect large quantum corrections 
to the hhh coupling constant. We note that these fermion loop contributions are always 
negative to the SM prediction. In FIG. [TBI effects of the chiral fourth generation fermions 
on the hhh coupling constant are shown for m^ = 120 GeV (left) and for m^ = 210 GeV 
(right). Hereafter, the mass differences are fixed to be m t / — my = 55 GeV with m t / = my 
and my = m v i. The masses of the fourth generation up-type quark are taken as three 
representative values, m t i = 256 GeV (dotted line), 300 GeV (dashed line), and 400 GeV 
(long-dashed line). The hhh coupling constant is changed significantly depending on the 
energy of the off-shell Higgs boson. In the low energy limit, a huge quantum correction to 
the hhh coupling constant can be more than 100%, which can easily overwhelm the SM 
contribution and change the sign of the total amplitude. We again note that from the 
vacuum stability condition it is highly disfavored to introduce too heavy fourth generation 
fermion unless the Higgs sector is extended. 

In FIG. [191 we show the invariant mass distribution of the differential cross section of 
gg — > hh in the SM with a complete set of fourth generation for m^ = 120 GeV (left) and for 
m/i = 210 GeV (right). The production cross sections of the chiral fourth generation model 
can be 10-100 times larger than that of the SM. The first peaks come not only from the 
large one-loop correction to the hhh vertex but also from the enhancement of the ggh vertex 
due to the fourth generation quarks. The threshold enhancements of on-shell fermion-pair 
production can also be seen around Mhh — 2m t and Mhh — 2m f which are smeared by the 
distribution of high energy gluons. 

In FIG. [201 we show the cross section of the double-Higgs-strahlung process e + e~ — > hhZ 
as a function of the collision energy. The cross section can be reduced by the suppression 
of the hhh coupling constant due to the fourth generation fermions. Unlike the case of 
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FIG. 19: The invariant mass distribution of gg — > hh process at the LHC with \J~s = 14 TeV for 
nih = 120 GeV (left) and nih = 210 GeV (right) in the chiral fourth generation model. Three 
representative values of t' mass are chosen as 256 GeV (dotted line), 300 GeV (dashed line), and 
400 GeV (long-dashed line). The SM prediction is also shown by a solid curved line for comparison. 
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FIG. 20: The cross section of e + e — > hhZ process as a function of y/s for mh = 120 GeV (left) 
and nih = 210 GeV (right) in the chiral fourth generation model. 



Ak approximation the deviation of the cross section from the SM depends on the collision 
energy, because the quantum corrections to the hhh coupling constant are the function of 
the energy for the off-shell Higgs boson. 

In FIG. EU the cross section of the Higgs pair production at a photon collider is given 
as a function of ■y/s^ with nih = 120 GeV (left) and nih = 210 GeV (right). All the fourth 
generation fermions contribute to the both triangular and the box diagrams (see FIG. [5]) 
due to the large Yukawa coupling constant, which can enhance the cross section significantly 
by a factor of 10 for the wide range of ^fs^ e . The threshold effects of the on-shell heavy 
fermions can be found soon above the thresholds ^J~s^ e ~ 2m t and ^/s^ ~ 2m f because 
the photon luminosity has a peak around 0.8^/5^. At the PLC, we can have larger cross 
sections even for relatively heavy Higgs bosons. 

In FIG. [221, we show the cross section for e + e~ — > hhvv with the one-loop corrected 
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FIG. 21: The cross section of 7(+)7(+) — > hh process at the photon collider option at the ILC as 
a function of the e~e~ collision energy for rrih = 120 GeV (left) and rrih = 210 GeV (right) in the 
chiral fourth generation model. 
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FIG. 22: The cross sections of e + e — > hhuv process at the ILC as a function of collision energy 
y/s for rrih = 120 GeV (left) and rrih = 210 GeV (right) in the chiral fourth generation model. 



hhh coupling constant. For rrih — 120 GeV (left), the production rate becomes significantly 
large compared to the SM rate. This enhancement mainly comes from the large quantum 
corrections to the hhh coupling constant in the smaller Mhh region. For rrih = 210 GeV, 
corrections to the hhh coupling constant are relatively small. However, for larger Mhh values 
the one- loop correction (AT^/pSMj g Qes k^r, ^ Q negative, and its effect rapidly becomes 
important for Mhh ^ 1500 GeV, which makes the cross section larger. 



D. Vectorlike quarks 



Various types of vectorlike fermions have also been discussed in the literature. They can 
appear in extra-dimension models with bulk fermions 69], in little Higgs models 70[ and 
in the top seesaw model 



71 



As a representative case of these models, we adopt a pair 
of vectorlike up-type quarks, T 0L and T 0R , which transform as (3, 1)2/3 under the gauge 
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symmetry. 

The Lagrangian relevant to the mass of the SM top-quark and the vectorlike up-type 
quark can be written as 



£-ass = _ yt Qq ^ <5>_ Yt Q T 0R $ - M T T 0L T 0R + H.c, 



(17) 



where we have dropped the terms proportional to T 0L t 0R which are absorbed by redefinitions 
of t 0R and T 0R without loss of generality. Since the t -T mixing term is allowed by the 
symmetry, to an d To are no longer mass eigenstates. The mass eigenstates t and T are 
determined by diagonalization of the mass matrix, 



yt v Y T v 
M T 



U 



t i m t 



u 



rriT 



R, 



(18) 



where 



x 



T 



x 





where X = L,R. 



(19) 



The direct search for the vectorlike quarks has been performed 64J. Their main pro- 
duction modes at hadron colliders are gg —> TT, so that the lower bound of the up-type 
vectorlike fermion is basically the same as the fourth generation up-type quark, m T > 256 
GeV. 

The vectorlike quarks are also severely constrained by electroweak precision data 72 . 
For U = 0, the experimental values for oblique parameters are S = —0.04 ± 0.09 and 
T = 0.02 ± 0.09 where rrih = 117 GeV is assumed. The contributions to the S parameter 
due to the vectorlike top-quark (T) are less than 0.004 for Mt > 1000 GeV, which can 
be neglected. The to ~ T Yukawa coupling constant is set to be Yt = 1 throughout this 
analysis. On the other hand, the T parameter is rather sensitive to the model parameters. 
The one (two) sigma bound on the lowest value of Mt is 1100 (1700) GeV for = 120 
GeV. For m h = 160 GeV, these constraints are slightly milder, M T > 1100(1500) GeV at 
1ct(2o") confidence level. 

The one-loop correction to the hhh coupling constant due to the vectorlike top-quark is 
evaluated as 



"pVec 
1 hhh 
pSM 
1 hhh 



N c m T 



37r 2 v 2 mf i ( 



3m 2 t 



m 2 /m T 



2 

mi, 



rrin 




where y^ s = c L [c R yt — s r Yt)- Although there is a m T enhancement factor, the correction 
to the hhh coupling constant can not be large. This is because a factor (l — y^v / \\/2mt)) 
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FIG. 23: The rates for one- loop contributions from the vectorlike top-quark T to the hhh coupling 
constant for m/, = 120 GeV (left) and for m/, = 160 GeV (right). The to _ ^o Yukawa coupling is 
taken to be Yj< = 1, and the gauge invariant mass parameter Mt is chosen as 1100 GeV (dotted 
line), 1500 GeV (dashed line) and 1700 GeV (long-dashed line), respectively. 
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FIG. 24: The invariant mass distribution of gg hh process at the LHC with -^/s = 14 TeV for 
irih = 120 GeV (left) and rrih = 160 GeV (right) in the SM with the vectorlike top-quark. 



is approximately expressed as (Ytv / (v^mr)) for large m T , so that the correction to the 
hhh coupling constant decouples as \/m\. In FIG. [231 the effects on the hhh coupling 
constant due to the vectorlike top-quark are shown. The stringent experimental bounds from 
the electroweak precision data impose that mass of T particle is heavy, m T > 1100 GeV. 
Therefore, there can be no significant nondecoupling effect on the hhh coupling constant. 

In FIG. rjH we show the invariant mass distribution for the cross section of gg —> hh 
at the LHC with yfs = 14 TeV in the model with the vectorlike top-quark. The vectorlike 
top-quarks give new contribution to the both triangular and box diagrams. However, the 
deviations of the cross section from the SM value are at most 5%, which can not be large 
because of their decoupling nature. In FIG. [251 the cross sections for the double-Higgs- 
strahlung process are shown in the model with the vectorlike top-quark. The effects of the 
vectorlike fermions only appear in the hhh coupling constant. The impact of the vectorlike 
top-quarks is quite small also in this process. In FIG. [261 we show the cross section of 
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FIG. 25: The cross section of e + e — > hhZ process as a function of -y/s for 
and mh = 160 GeV (right) in the SM with the vectorlike top-quark. 
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FIG. 26: The cross section of t(+)t(+) — > hh process at the photon collider option at the ILC as 
a function of the e~e~~ collision energy for = 120 GeV (left) and = 160 GeV (right) in the 
SM with the vectorlike top-quark. 



77 — > hh process in the model with vectorlike top-quark. Similarly to the gluon fusion 
process gg — > hh, the new physics effects are rather small. In FIG. [27J the cross section for 
e + e~ — > hhvv is shown as a function of e + e~ energy. Since the effect on the hhh coupling 
constant are small in this model, the deviation of cross section is very tiny. We note that the 
deviations of cross sections for gg — > hh and 77 — > hh are larger than those for e + e~ — > hhZ 
and e + e~ — > hhvv due to the flavor changing Yukawa interaction between t and T from the 
one-loop box diagrams. 



IV. SUMMARY AND DISCUSSIONS 



We have studied the double Higgs boson production processes gg — > hh, e + e~ — > hhZ, 
77 — > hh and e + e~ — > hhvv in various new physics models. These processes include diagrams 
that contain the hhh coupling constant, so that they can be used to obtain information of 
the Higgs potential. SM cross sections for these processes are shown at the leading order in 
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FIG. 27: The cross sections of e + e — > hhuv process at the ILC as a function of collision energy 
y^s for m/j = 120 GeV (left) and = 160 GeV (right) in the SM with vectorlike top-quark. 



m h [GeV] 


4& h » 








120 


21 


0.09-0.2 


0.14 


0.09-0.9 


160 


12 


0-0.05 


0.11 


0.02-0.6 


185 


8.0 


0-0.01 


0.07 


0.02-0.5 


210 


5.2 





« 


0.01-0.4 



TABLE II: The total cross sections of Higgs boson pair production for = 120, 160, 185 and 210 
GeV in the SM are listed. The proton-proton collision energy is taken as \fs = 14 TeV for the 
gluon fusion mechanism. For the double-Higgs-strahlung - v /s^T is varied from 400 GeV to 500 GeV, 
while for the W boson fusion ^/s ee is varied from 1 TeV to 3 TeV. For photon-photon collisions 
energy is optimized to obtain the largest cross sections in a range ^Js ee > 500 GeV. 

Table Hll for several values of nth- In order to see the impact of a deviation in the hhh coupling 
constant on these double Higgs boson production processes, we have at first evaluated their 
cross sections in the SM but assuming the constant deviation in the hhh coupling constant 
by the factor of (1 + An). The results are summarized as follows: 

• The gg — > hh process is the one-loop induced process, where contributions from the 
triangle top-loop diagrams with the hhh coupling constant and the box-type top- 
loop diagrams are destructive. Therefore, a negative (positive) deviation in the hhh 
coupling constant can make the cross section larger (smaller). For example, when Ak = 
— 1, namely the case without the hhh coupling, the cross section can be approximately 
doubled (tripled) as compared to the SM value with An = for m h = 120 GeV (160 
GeV). 

• For the process of e + e~ — > hhZ, on the other hand, the contribution of the tree-level 
diagram with the hhh coupling constant and that of the other tree-level diagrams are 
constructive, so that the cross section is enhanced by the positive deviation in the 
hhh coupling constant. Because of an s-channel process, as seen in Table Hit the cross 
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section becomes rapidly smaller for a larger mass of the Higgs boson, so that this 
process may only be useful for a light Higgs bosons like < 140 GeV for yfs ~ 500 
GeV. 

• The one-loop induced process 77 — > hh can play a complementary role to gg — > hh 
and e + e~ — > hhZ. For -y/s^ < 500 GeV, the sensitivity to the deviation in the hhh 
coupling constant is significant for both = 120 GeV and 160 GeV. Direction of 
interference of the diagram with the hhh coupling constant and the other diagrams 
is however opposite; i.e., destructive and constructive for m h = 120 GeV and 160 
GeV, respectively. These characteristic behaviors of 77 — > hh can be complementary 
to gg — > hh and e + e~ — > hhZ in the measurement of the hhh coupling constant. 
Furthermore, sensitivity to the deviation in the hhh coupling constant can be better 
by using this process even when the collision energy of the linear collider is limited to 
be relatively low (^/~sTe ^ 500 GeV). 

• The W boson fusion process e + e~ — > hhvv can be useful for measuring the hhh cou- 
pling constant at an energy upgrade of the ILC or the CLIC with ^fs^ e = 1-3 TeV, 
because the production cross section is monotonically increasing with the collision en- 
ergy due to the t-channel colinear effect. The diagram with the hhh coupling constant 
has the opposite sign with the other diagrams. 

We have evaluated cross sections for these processes in various new physics models such 
as the THDM, the model with scalar leptoquarks, the model with chiral fourth generation 
quarks, and the model with a vectorlike quark. In these models, apart from the deviated 
hhh coupling constant, additional one- loop diagrams can contribute to the cross sections 
especially in the loop induced processes gg —> hh and 77 — > hh. Only one-loop diagrams 
of colored particles contribute to the former process, while those of all the charged particles 
do to the latter one. In Table IHIt the results for possible deviations in cross sections for 
these processes are summarized with the deviation in the hhh coupling constant for several 
typical values of mh and the collision energies in each model. We summarize the results for 
each model below in order. 

In the THDM, one-loop corrections of additional scalar bosons to the hhh coupling con- 
stant can give +100% deviations in the SM-like limit where we set sin(/3 — a) = 1 with 
M 2 ~ 0. This positive large quantum correction to the hhh coupling constant is common 
in the extended Higgs sectors with nondecoupling property where the mass of the scalar 
bosons comes mainly from the VEV. The effect of the large deviation in the hhh coupling 
constant can be well described by the analysis with the constant shift of the hhh coupling 
constant by the factor of (1 + An). A qualitative difference can be seen in the one- loop 
induced 77 — > hh process, where one-loop diagrams of charged Higgs bosons can change to 
the cross section. 

In the scalar leptoquark models, the correction to the hhh vertex is positive because 
of additional bosonic loop contributions, as in the THDM. The magnitude can be larger 
than +100% via the nondecoupling effect of scalar leptoquarks in the loop when Mjg ~ 0. 
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Model 


m h [GeV] 


N 1 iJ \ 

T^IN F T^b 1V1 

hhh hhh 
pSM 
hhh. 




A„,e+e ~>hhZ 


^ r NP 


A„e+e —>hhvv 


THDM 


120 


+120% 


-50% 


+(80-70)% 


+50% 


-(80-50)% 


THDM 


160 


+70% 


-50% 


+ (60-50)% 


+110% 


-(80-50)% 


LQ(g = 1/3,4/3) 


120 


+150% 


-40% 


+(110-100)% 


+130%, +100% 


-(70-60)% 


LQ(Q = 1/3,4/3) 


185 


+60% 


-30% 


+50% 


+150%, +150% 


-(80-50)% 


Ch4 


120 


-590% 


+7800% 


-(30-20)% 


+3100% 


+(260-110)% 


Ch4 


210 


-140% 


+2200% 






+(970-210)% 


Vec 


120 


-4% 


-10% 


-2% 


-10% 


+(5-l)% 


Vec 


160 


-2% 


-5% 


-1% 


-10% 


+(3-0)% 



pNP (4 m 2 m 2 m 2\ 

TABLE III: Possible quantum corrections to the hhh coupling constant, r sM ( d 2' — t — T\ ~ 1) an d 

hhh^ m h' m h> m h' 

deviations of cross sections Ar NP = (ctnp — csmV^sm are listed. The proton-proton collision energy 
is taken as y/s = 14 TeV for the gluon fusion mechanism. For the double-Higgs-strahlung y/s ee is 
varied from 400 GeV to 500 GeV, while for the W boson fusion y/s ee is varied from 1 TeV to 3 
TeV. For photon-photon collisions y/s ee is optimized to obtain the largest cross sections. Model 
parameters are chosen as m$ = 400 GeV and M 2 = for THDM, = 400 GeV and Mjn = 
for scalar leptoquark models, m? = 400 GeV and my = 345 GeV for the fourth generation model, 
and Yt = 1,Mt = 1100 GeV for vectorlike top-quark model, respectively. 



A qualitative difference from the THDM case is that the scalar leptoquarks are colored, 
which can contribute to the gg — > hh through the one particle irreducible one-loop diagram. 
However, it turns out that the top-quark one-loop contribution is much larger than the 
leptoquark-loop contribution, so that the SM result with An correction is a good approxi- 
mation. Therefore, as expected in the analysis by using the An, positive deviations in the 
hhh coupling constant make the cross sections smaller. It amounts to minus 40(30)% for 
rrih = 120(185) GeV assuming the other parameters as = 400 GeV and Mj^ = 0. The 
production rates for e + e~ — > hhZ can be enhanced by +100(+50)% for rrih = 120(185) GeV 
due to the constructive interference in the contribution from the diagram with a positively 
deviated hhh coupling constant and the other diagrams. On the contrary, cross sections for 
e + e~ —> hhvv become smaller due to the destructive interference. The production rates for 
77 — > hh depend on electric charges of leptoquarks. For the scalar leptoquark with Q = 4/3, 
the cross section can be enhanced by the threshold effect at ^/s 77 ~ 2m^, while for those 
with Q — 1/3 such effects are smeared. For yfs^ > 500 GeV, the cross section for the scalar 
leptoquarks with Q = 4/3 can enhance more than several times +100%. 

In the model with chiral fourth generation quarks, the hhh coupling constant can be 
changed by more than a few times —100% due to the nondecoupling loop effect of additional 
heavy chiral fermions. These huge corrections can be possible under the constraint from 
the data for precision measurements and the direct search results. However, such a large 
fermionic loop contribution can make Higgs potential unstable, so that a heavier Higgs 
boson is required than the allowed value in the SM to recover the stability of vacuum. 
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A light Higgs boson can be allowed by extending the Higgs sector with additional scalar 
doublets. The cross section for gg — > hh is drastically enhanced from the SM prediction 
by a factor of 10-100, because new colored particles contribute to the additional one-loop 
diagrams at leading order and the fourth generation fermions enhance both the one-loop 
induced vertex ggh and the one-loop corrected hhh coupling constant. For a reference 
point of m t i = rri£ = 400 GeV and my = m u i = 345 GeV, the cross sections become 7800 
(2200)% for m h = 120 (210) GeV. Consequently, the process gg — > hh with the decay mode 
of h — > WW/ZZ can be promising to measure the hhh coupling constant for > 210 
GeV. For e + e~ — > hhZ process, the effect of fourth generation fermions only appear in the 
one-loop corrected hhh coupling constant. Since this correction is negative because of the 
fermionic loop contribution, the cross section is suppressed. Numerically, it is reduced by 
30-20% for the above reference point with m h = 120 GeV and 400 GeV < Js^ e < 500 
GeV. On the contrary, in e + e~ — > hhvv process, the effect of fourth generation fermions 
only appears in the hhh coupling constant through large one-loop corrections, which makes 
the cross section huge as compared to the SM value. Cross sections for 77 — > hh at photon 
colliders can also be modified by extra fermion loops similarly to those for gg — > hh. With 
an optimized value of ^fs^ e < 500 GeV, the cross section can be enhanced by 3100% for 
the reference point with = 120 GeV. Since the cross sections are drastically enhanced 
for rrih = 120 GeV and hadronic decay modes can be measurable at the PLC, the process 
77 — > hh would be a promising process to probe the hhh coupling constant. 

In the vectorlike quark model, a nondecoupling limit (My ~ 0) cannot be taken due to the 
severe experimental constraints. Therefore, there are no large one-loop effect from a vector 
top-quark on the hhh coupling constant. The cross sections for gg — > hh and 77 —> hh 
can be deviated slightly from SM prediction due to the flavor changing Yukawa interaction 
of vectorlike quarks. However, deviations of Higgs boson pair production cross sections are 
rather small, so that huge luminosity would be required for measuring the deviation of the 
hhh coupling constant. 

Measuring four kinds of double Higgs boson production processes at different future col- 
lider experiments is useful to discriminate whether new physics particles in the loop are 
fermions or bosons and also whether they are colored or not. Higgs boson pair production 
processes e + e~ —> hhZ and e + e~ — > hhvv at lepton colliders are tree-level processes, which 
enable us to extract information of the hhh coupling constant. On the other hand, measure- 
ments of effective vertices gghh and 77/1/1 in loop induced processes gg — > hh at the LHC 
and 77 — > hh at the PLC can provide information of colored and electrically charged parti- 
cles in loop diagrams. Effective vertices ggh and 77/1 can be determined in the single Higgs 
boson production processes gg — > h and 77 — > h as well. Combining these measurements, 
we would be able to disentangle new physics effects in the hhh coupling and the effective 
vertices. 

We have considered the Higgs boson pair production processes, gg — > hh, e + e~ —> hhZ, 
e + e~ — > hhvv and 77 — > hh as a probe of the hhh coupling constant. The measurement 
of the hhh coupling constant is particularly important to understand the mechanism of the 
electroweak symmetry breaking. The hhh coupling constant can receive quite large quantum 
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corrections from new physics particles as a nondecoupling effect, which can be an order of 
more than 100%. Deviations of the hhh coupling constant can give different effects on 
these processes which can largely modify production cross sections. Additional particles in 
new physics model can also significantly affect the gg — > hh and the 77 — > hh processes 
according to their color and electric charges. We have found that these four Higgs boson 
pair production processes at different colliders can play complementary roles in exploring 
new physics through the Higgs sector. 
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Appendix A: The one-loop corrections to the triple Higgs boson coupling 



The relatively large one-loop correction to the hhh coupling constant in the SM has been 
calculated as 
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where B and C are the loop functions, which are defined in Ref. 73]. For the chiral fourth 
generation model, all the extra fermions further contribute to the hhh coupling constant. 
Then the sum should be replaced by all heavy fermions in the preceding formula. 
In the leptoquark model, the quantum effect is given by 
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In the THDM with the SM-like limit sin(/3 — a) = 1, the one-loop corrections to the hhh 
coupling constant are calculated as 
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In the vectorlike top-quark model, the one-loop corrected hhh coupling constant is given 
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+ (yf^yf,m t ^m T )^, (A7) 
where 

yf=c L (c R y t -s R Y T ), (A8) 

e t = s L( c R.yt - s R Y T ), (A9) 

y e T=s L (s R y t + c R Y T ), (A10) 

tT= C L( S Ryt + C R Y T ). (All) 



Appendix B: The loop integrals for gg — > hh and 77 — > hh 



The SM contributions to the loop functions in gg — > hh and 77 — > amplitudes can 
be found in Refs. 
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291 ] . It can be generalized straightforwardly for the chiral fourth 



generation model. 

In leptoquark models, in addition to the SM fermions there are contributions from the 
colored scalar particle to gg — > hh amplitude as, 
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and the loop functions are abbreviated as 

y(1,2) 



C A[i,j,k] = C a( q , 0, s; mj,m k ), 

,(3,4) 
A[ij,fc] 



= C A (s, m^; m*, my, m fc ) 



(B5) 
(B6) 



33 



^(1,3) 
U A[i,j,k] 

^(2,3) 
^A[i,j,k] 
n (l,2,3) 
u A\ij,k,l\ 



D 



(2,1,3) 
A[i,j,k,l] 

n (l,3,2) 



C A (i, 0, m 2 h ; im, rrij, m k ), 

D A (0, 0, ml, m 2 h , s, u; m i; mj, m k , mi), 
D A (Q, 0, ml, m 2 h , s, i; m u m jy m k , mi), 
D A (0, m 2 h , 0, m 2 h , i, u; m u mj, m k , m t ). 



(B7) 
(B8) 
(B9) 
(BIO) 
(Bll) 



The index A denotes classes of loop functions, and we omitted obvious mass indices of 
intermediate particles in the Eqs. (lBlj) - (lB3j) . The leptoquarks also contribute to 77 — > hh 
process. These additional contributions to the amplitude are basically the same as those 
loop integrals for gg — > hh amplitude except for electric charges and a color factor; = 
N+Q\F* and H+~ = N*Q\G*. 

The SM contribution can be modified in the vectorlike quark model because of t and 
T mixing. The vectorlike top-quark T also gives additional contributions to gg — > hh and 
77 — > hh amplitudes as 
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